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ABSTRACT 


Heat  exchangers  with  one  shell  pass  and  n  tube  passes 
are  often  referred  to  as  1-n  exchangers.  The  heat  transfer 
literature  contains  many  references  to  studies  of  1-n 
exchangers  when  n  is  even  but  apparently  little  work  has 
been  done  with  respect  to  the  1-n  exchanger  when  n  is  odd. 
This  thesis  greatly  expands  the  theoretical  study  of  1-n 
exchangers  with  n  being  odd.  While  a  completely  closed  form 
solution  was  found  to  be  unfeasible,  a  polynomial  approxi¬ 
mation  has  been  developed  that  yields  the  effectiveness  ( e) 
of  the  two  possible  arrangements  of  the  1-3  exchanger  as  a 
function  of  the  capacity  rate  ratio  (R)  and  the  number  of 
transfer  units  (Ntu).  It  is  also  shown  that  the  effective¬ 
ness  of  the  arrangement  with  two  counterflow  and  one 
parallel  flow  tube  side  passes  exceeds  that  of  some  of  the 
1-n  exchangers  with  n  even. 
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NOMENCLATURE 


English  Letter  Symbols 

A  *  Exchanger  heat-transfer  surface,  sq  m 

Am  *  Coefficient  of  the  rath  value  dimensionless 

Ao  53  Coefficient,  dimensionless 

Ai  *  1st  order  coefficient  to  be  multiplied  by  Ntu, 
dimensionless 

A2  *  2nd  order  coefficient  to  be  multiplied  by  Ntu2, 
dimensionless 

A3  *  3rd  order  coefficient  to  be  multiplied  by  N-tu3 , 
dimensionless 

A4  *  4th  order  coefficient  to  be  multiplied  by  N-^4 , 
dimensionless 

A5  *  5th  order  coefficient  to  be  multiplied  by  Ntu5, 
dimensionless 

a  =  Exchanger  heat-transfer  surface,  sq  m/m 

C  =  Capacity  rate,  W*/K.  Also  designates  dimensionless 
arbitrary  constant 

Cpc  ■  Specific  heat  at  constant  pressure  of  cold  fluid, 
J/kg-K 

Cph  *  Specific  heat  at  constant  pressure  of  hot  fluid, 

J /kg°K 

D  *  Empirical  value  of  effectiveness  (computer  generated), 
dimensionless 

e  =  Error.  Also  used  as  the  exponential  function 

F  =  Logarithmic  mean  temperature  difference  correction 
factor,  dimensionless 

L  *  Exchanger  length,  m 

N  *  Number  of  effectiveness  empirical  data  points  used  to 
determine  a  curve  for  R,  dimensionless 


n  ■  Number  of  tube  passes,  dimensionless.  Also  number  of 
equations,  dimensionless 

nc  -A  related  Ntu  per  unit  length  hot  side,  m”1 
nh  “A  related  Ntu  per  unit  length  cold  side,  m"1 
Ntu  “  Number  of  transfer  units,  dimensionless 
P  *  Temperature  group,  dimensionless 
q  *  Total  rate  of  heat  transfer,  W 
Qmax=  Maximum  total  rate  of  heat  transfer,  W 
R  *  Capacity  rate  ratio,  dimensionless 
S  ■  Temperature  group,  dimensionless 

Sr  =  Sum  of  the  squares  of  the  residuals,  dimensionless 

T  *  Hot  fluid  temperature,  *C 

Tpi  ■  Particular  integral,  dimensionless 

Ti  =  Hot  fluid  temperature  in,  *C 

T2  *■  Hot  fluid  temperature  out ,  *C 

ti  *  Cold  fluid  temperature  in,  °C 

t2  ■  Cold  fluid  temperature  out,  °C 

ta  *  Cold  fluid  temperature  1st  pass,  *C 

tab  *  Cold  fluid  temperature  between  1st  and  2nd  passes 

tb  38  Cold  fluid  temperature  2nd  pass,  °C 

tbc  “  Cold  fluid  temperature  between  2nd  and  3rd  passes 

tc  *  Cold  fluid  temperature  3rd  pass,  *C 

U  »  Overall  heat  transfer  coefficient,  W/m2  -  #C 

W  =  Mass  flow,  kg/sec.  Also  the  product  of  u  and  L, 

dimensionless 

x  *  length  coordinate,  m.  Also  used  to  represent 
a  constant  value  in  a  sequence 


y  -  Sum  of  mth  degree  polynomial,  defined  by  eq.  (53), 
dimensionless 

Z  ■  A  product  of  z  and  L,  dimensionless 
z  -  A  related  N^u  per  unit  length,  hot  side,  1/m 

Greek  Letter  Symbols 

a  *  Root  of  auxiliary  differential  equation,  1/m 

e  ■  Exchanger  effectiveness,  dimensionless 

X  -  Combination  of  variables  defined  by  equation  (11), 
dimensionless 

CD  *  A  related  Ntu  per  unit  length,  hot  side,  1/m 

4>  *  A  combination  of  terms  defined  by  eq.  (38), 

dimensionless 

Z  *  Summation,  dimensionless 

02  s  Variance,  dimensionless 

em  *  Mean  temperature  difference  for  exchanger,  *C 
3  ■  Indicates  partial  derivative,  dimensionless 

Subscripts 

c  *  Cold  fluid 

h  =  Hot  fluid 

i,j,k  *  Values  in  a  sequence 

m  =  Degree  or  order,  an  exponent 

1  *  inlet 

2  »  outlet 

Special  Symbols 

[A]  «  An  m  x  n  matrix,  symmetric 
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I .  INTRODUCTION 


A.  BACKGROUND 

When  analyzing  the  standard  counterflow  heat  exchanger, 
it  becomes  apparent  that ,  from  a  practical  standpoint ,  it  is 
often  difficult  to  obtain  a  high  velocity  for  one  of  the 
fluids  when  this  fluid  is  constrained  to  flow  through  all  of 
the  tubes  in  a  single  pass.  This  leads  to  a  possibility  of 
a  low  overall  heat  transfer  coefficient  which  cancels  the 
advantage  of  the  high  logarithmic  mean  temperature 
difference  which  is  obtainable  in  true  counterflow. 

The  quest  for  flow  arrangements  for  Increased  heat 
recovery  has  led  to  arrangements  that  yield  increased  tube- 
side  velocities  and  higher  overall  heat  transfer 
coefficients  even  at  the  expense  of  a  departure  from  the 
ideal  true  counterflow  arrangement.  Thus,  the  design  may  be 
modified  so  that  the  tube  side  fluid  is  carried  through 
fractions  of  the  tubes  consecutively. 

Heat  exchangers  of  this  type  with  one  shell  pass  and  n 
tube  passes  are  often  referred  to  as  1-n  exchangers.  These 
exchangers,  such  as  the  one  shell  pass,  two-tube  pass  (1-2) 
parallel-counterflow  exchanger  (see  Figure  1.1),  are 
configured  such  that  all  of  the  tube  side  fluid  flows 
through  the  two  halves  of  the  the  tubes  successively.  A 
single  channel  is  employed  with  a  partition  to  permit  the 


Figure  1.1  1-2  Parallel-Counterflow  Exchanger 

entry  and  exit  of  the  tube  side  fluid  from  the  same  channel. 
Note  the  baffles  used  to  induce  turbulence  causing  the 
liquid  to  flow  through  the  shell  at  right  angles  to  the  axes 
of  the  tubes  thus  helping  to  create  a  higher  shell  side 
velocity  with  higher  shell  side  heat-transfer  coefficients. 

To  date,  much  work  has  been  done  on  finding  the  true 
logarithmic  mean  temperature  difference  for  heat  exchangers 
with  an  even  number  of  tube  passes.  Little  work,  however, 
has  been  done  with  regard  to  heat  exchangers  having  an  odd 
number  of  passes.  This  is  primarily  due  to  the  method 
employed  in  deriving  an  analytical  solution  to  measure  the 
overall  effectiveness  of  a  heat  exchanger. 


Exhangers  with  an  even  number  of  tube  passes  often 
present  a  configuration  problem  especially  in  a  marine  appli 
cation  where  the  inlet  and  outlet  of  the  cooling  fluid  must 


be  one  the  same  side  of  the  exchanger  header  (see  Figure 
1.1).  This  particular  problem  could  be  alleviated  by  going 
to  heat  exchangers  with  an  odd  number  of  tube  passes. 
Currently,  precise  mathematical  expressions  for  the 
effectiveness  of  the  1-3  exchangers  do  not  exist.  Hence,  a 
theoretical  examination  is  reported  on  here  which  considers 
the  effectiveness  of  the  1-3  parallel-counter  flow  exchanger 
which  is  shown  in  Figure  1.2. 

Before  doing  this  it  is  important  to  mention  the  work 
that  lead  to  the  effectiveness  method  and  the  development  of 
work  on  heat  exchangers  with  an  odd  number  of  tube  passes. 


Figure  1.2  1-3  One  Shell  Pass  Three  Tube  Pass  Exchanger 


This  will  be  considered  in  Section  II.  Kern  [Ref.  1: 
pp.  224-226],  makes  the  interesting  point  that  the  optimum 
exchanger  requires  an  exchanger  capable  of  providing  the 
optimum  fluid-flow  velocities  on  the  shell  as  well  as  the 
tube  sides.  This  might  frequently  entail  the  use  of  an  odd 
number  of  tube  passes  or  an  odd  tube  length. 


B.  WHY  EFFECTIVENESS  AS  A  FUNCTION  OF  Ntu 

It  is  also  noted  that  Kays  and  London  [Ref.  2:  pp. 

24-29]  indicated  that  the  effectiveness  as  a  function  of 
Ntu  (e  "  Ntu)  method  is  the  favored  approach  for  evaluating 
a  heat  exchanger's  performance  because: 

1.  The  effectiveness  value  stands  alone  as  a  dependent 
variable  and  should  not  appear  directly  in  the 
abscissa  and  indirectly  in  the  ordinate  of  a  graphical 
display. 

2.  The  log-mean  difference  equation  misleadingly 
simplifies  the  notion  of  what  is  involved  in  heat 
exchanger  design  theory,  since  the  implication  is  that 
only  a  rate  equation  is  required. 

3.  The  e  -  Ntu  approach  simplifies  the  algebra  involved 
in  predicting  the  performance  of  complex  flow 
arrangements. 

4.  The  more  meaningful  arguments  are  related  to  ease  of 
use  in  design  work.  Two  prime  examples  of  these  are: 

a)  Given  the  overal  heat  transfer  coefficient,  U,  the 
two  fluid  capacity  rates,  Cc  and  C^,  and  the 
terminal  temperatures,  determine  the  required 
surface  area,  A. 

b)  Given  A,  U,  Cc,  Ch  and  the  inlet  temperatures  of 
both  streams,  determine  the  outlet  temperatures. 


II.  THE  DEVELOPMENT  OF  THE  EFFECTIVENESS  METHOD 

A.  LITERATURE  SURVEY 

It  was  Nagle  [Ref.  3:  pp.  604-609],  in  1931,  who 
credited  Davis  [Ref.  4]  with  a  simplified  method  for 
computing  actual  temperature  differences  between  two  heat¬ 
interchanging  streams  which  depart  from  true  counter  or 
concurrent  (parallel)  flow.  This  is  now  the  familiar  "F 
factor"  method  which  expresses  the  actual  mean  temperature 
difference  0m  in  q  =  UAem  as  a  fraction  F  of  the  counterflow 
logarithmic  mean  temperature  difference,  LTMD,  e^c  via  0^  = 

F  ®mc  • 

The  example  of  initial  interest  was  the  1-2  exchanger 
with  a  single  shell  pass  and  two  continuous  tube  passes  in 
counter  and  concurrent  flow  with  it.  The  method  involved 
derivation  of  the  actual  temperature  difference  for  the  flow 
pattern  and  formed  the  ratio  F  *  9m/9mc*  This  familiar  LMTD 
correction  factor  was  plotted  conveniently  as  functions  of 
the  effectiveness,  e,  and  the  capacity  rate  ratio  P  with  R 
as  a  parameter.  These  mean  temperature  difference 
correction  charts  are  available  for  many  flow  arrangements 
fRef.  1:  pp.  829-833  and  Ref.  5].  The  effectiveness,  e 
(often  called  P  or  S),  is  always  the  cold  fluid 
effectiveness  and  R  is  always  the  capacity  rate  ratio  of 


cold  fluid  to  hot  fluid. 


Nagle  detailed  assumptions  and  derivations  for  the  1-2, 
1-4  and  1-6  exchangers.  The  P  factors  were  obtained  by 
Nagle  through  graphical  integration  and  were  accomoanied  by 
the  comment  that  F  factors  for  the  1-2  exchanger  could  be 
applied  with  negligible  error  to  1-4  and  1-6  exchangers. 
Underwood  [Ref.  6:  pp.  145-148]  rederived  the  equations  of 
Nagle  for  1-2  and  1-4  exchangers  to  eliminate  the  need  for 
obtaining  F  factors  by  graphical  integration. 

Bowman  [Ref.  7:  pp.  541-544  ]  pointed  out  that  for  a 
very  large  or  infinite  number  of  tube  passes,  the  F  factor 
approached,  as  a  limit,  its  value  in  crossflow  with  both 
fluids  completely  mixed.  It  was  further  stated  that  even  at 
the  limit,  the  F  factors  were  only  1  to  2  percent  lower  than 
those  for  the  1-2  exchanger.  A  previous  paper  by  Kraus  and 
Kern  [Ref.  8]  did  not  confirm  the  generalization  that  1-n 
exchangers  differed  only  negligibly  from  the  1-2  exchanger 
although  this  lack  of  confirmation  was  obtained  on  an 
e  *  f(R,  N-tu )  basis.  Moreover,  the  Kraus-Kern  work  does  not 
confirm  the  generalizations  on  an  F  *  f(R,  Ntu>  e)  basis. 

From  the  standpoint  of  usefulness  and  good  accuracy,  it 
is  essential  that  F  factors,  if  they  are  to  be  used  in 
nreference  to  e  =  f(R,  N^u»  flow  arrangement),  he  obtained 
with  precision.  Plots  of  F  *  f(R,  e,  *  P  or  S)  [Ref.  1: 
pp.  829-833  and  Ref.  5]  show  that  the  curves  for  particular 
values  of  R  approach  infinite  slope  as  F  decreases.  While 
this  can  be  partially  alleviated  by  restricting  R  <  1.0  (a 
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constraint  used  in  the  e  *  f(R,  N-tu*  flow  arrangement) 
approach),  it  is  seen  that  small  errors  in  the  interpolation 
for  R  or  e  »  P  or  S  can  result  in  large  fluctuations  in  the 
value  of  F. 

In  a  comprehensive  paper.  Bowman,  Mueller  and  Nagle 
[Ref.  9:  pp.  283-294]  presented  graphs  of  F  factors  for 
shells  with  one  through  six  shell  passes  and  numbers  of 
continuous  tube  passes  respectively  double  the  number  of 
shell  passes.  In  view  of  the  earlier  references  to  Nagle 
and  Bowman,  it  should  be  noted  that  F  factors  were  computed 
for  the  1-2  exchanger  in  [Ref.  9:  pp.  283-294]  using  the 
equations  of  Underwood  [Ref.  6:  pp.  145-158]. 

Ten  Broeck  [Ref.  10:  pp.  1041-10421  prepared  a  graph 
of  the  dimensionless  groups  now  known  as  e,  R  and  N^u  for 
the  1-2  exchanger.  Such  a  graph  had  the  added  versatility 
of  simplifying  the  calculation  of  performance  in  a  given 
exchanger  when  operating  at  conditions  different  for  those 
for  which  it  was  designed.  Kays  and  London  [Ref.  2: 
pp.  63-74]  prepared  similar  graphs  and  tables  of 
e  *  f(R,  Ntu*  flow  arrangement)  for  the  1-2  exchanger  and 
for  several  cases  of  crossflow  and  periodic  flow. 

The  foregoing  describes  the  early  history  of  the  search 
for  the  so-called  Logarithmic  Mean  Temperature  Difference 
Correction  Factor,  F,  with  regard  to  heat  exchangers  having 
an  even  number  of  tube  passes.  It  is  a  fact,  however,  that 
certain  space  economies  could  be  realized  from  exchangers 
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having  an  odd  number  of  tube  passes  so  that  the  tube  side 
fluid  could  enter  and  leave  the  exchanger  at  opposite  ends 
of  the  exchanger  (see  Figure  1.2). 

B.  FISCHER'S  WORK 

With  the  foregoing  in  mind,  an  extensive  search  has  been 
conducted  to  obtain  e  -  N-^u  data  for  the  so  called  "1-3"  and 
"1-5"  exchangers.  This  search  has  uncovered  a  single  work, 
that  of  Fischer  fRef.  11:  pp.  377-383),  which  summarizes 
the  historical  development  covered  here  and  contains  only  a 
small  section  on  the  1-3  exchanger.  This  work  by  Fischer 
develops  an  equation  for  true  mean  temperature  difference  of 
the  1-3  exchanger  and  casts  the  results  in  terms  of  F  rather 
than  e.  Moreover,  the  work  treats  only  the  case  where  the 
three  tube  passes  are  arranged  with  two  in  counterflow  and 
one  in  parallel  flow  (1-3:2C)  making  no  mention  of  the  one 
counterflow  and  two  parallel  flow  (1-3:2P)  case  (see  Figures 
2.1  and  2.2).  In  addition,  the  equation  developed  to  yield 
F  must  be  solved  using  a  trial  and  error  solution. 

The  present  work  is  aimed  at  continuing  the  Fischer 
investigation  for  several  reasons: 

1.  A  solution  is  needed  for  effectiveness,  e,  as  a 
function  of  capacity  rate  (R)  and  number  of  transfer 
units  (Ntu). 

2.  This  solution  should  he  in  a  closed  form  if  at 
all  possible  so  that  it  will  be  computationally 
efficient  and  useful  in  both  the  design  and  analysis 
frameworks. 
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3.  Because  valid  design  data  can  evolve  from  a 
polynomial  approximation.  The  search  should  not  be 
abandoned  just  because  a  closed  form  solution  does  not 
result  from  an  analytical  approach. 

4.  Data  is  needed  for  the  (1-3:2P)  two  parallel-one 
counterflow  configuration. 

5.  A  1-3  exchanger  in  a  marine  (shipboard)  application 
may  result  in  a  considerable  space  saving  over  its 
1-n  counterpart  with  n  even.  This  would  be  evident 
on  the  outside  of  the  exchanger  where  it  would  be 
immediately  noted  that  the  1-3  exchanger  has  tube 
side  inlet  and  outlet  at  opposite  ends  of  the 
exchanger. 

The  next  section  confirms  Fischer's  result  and  shows 
that  a  closed  form  solution  cannot  be  obtained  for  the 
effectiveness  of  the  1-3  exchanger.  Sections  IV  and  V 
demonstate  how,  through  numerical  analysis  assisted  by  a 
computer,  a  polynomial  solution  can  be  derived  that  will 
yield  the  effectiveness  to  engineering  accuracy. 


III.  AN  ATTEMPT  AT  A  CLOSED  FORM  SOLUTION 


A.  EFFECTIVENESS  AS  A  FUNCTION  OF  CAPACITY  RATES  AND 
EXCHANGER  SIZE 

This  section  deals  with  an  investigation  into  the 
effectiveness,  e,  of  a  one  shell  pass  and  three  tube  pass 
heat  exchanger,  whereby  e  compares  the  actual  heat  transfer 
rate  to  the  thermodynamically  limited,  maximum  possible  heat 
transfer  rate  as  would  be  realized  only  in  a  counter  flow 
heat  exchanger  of  infinite  transfer  area.  This  exchanger 
heat  transfer  effectiveness  is  given  by 

q  ^  ^h^hot^n"  Thot,out^  m  ^c^cold .out”  tcold,in^ 
qmax  Cmin^hot,in”  * cold, in ^  Cmin^hot,in~  cold, in) 

where  Cmin  is  the  smaller  of  the  C^  and  Cc  magnitudes. 

Thus,  e  possesses  the  significance  of  effectiveness  of  the 
heat  exchanger  from  a  thermodynamic  point  of  view,  with  the 
magnitude  of  the  effectiveness  completely  defining  the  heat 
transfer  performance.  In  general  we  express  e  *  f(Ntu,  R, 
and  flow  arrangement)  and  when  the  flow  arrangement  is  under¬ 
stood,  it  is  said  that  e  *  f(Ntu*  R)*  fRef.  2:  pp.  14-26). 

The  number  of  heat  transfer  units  Ntu  is  a  nondimensional 
expression  of  rhe  "heat  transfer  size"  of  the  exchanger. 

When  N^u  is  small  the  exchanger  effectiveness  is  low,  and 
when  N^u  is  large,  e  approaches  the  limit  imposed  by  the  flow 


arrangement  and  thermodynamic  conditions  asymptotically 
From  inspection  of  the  definition  of  N-tu 


it  is  clear  that  the  overall  conductance  and  transfer 
area  affect  the  costs  of  attaining  a  high  value  for  Ntu, 
ergo  high  e.  The  capacity  rate  ratio,  R,  as  defined  by 


is  simply  the  ratio  of  mass  flow  rate  times  specific  heat 
capacity  for  the  two  streams.  These  can  be  considered  as 
flow  stream  thermal-capacity  rates,  i.e.,  energy  storage 
rate  in  the  stream  per  unit  of  temperature  change.  [Ref.  2: 
pp.  14-26] 

The  attempt  taken  in  this  thesis  to  develop  a  closed 
form  solution  has  used  the  basic  fundamentals  of  heat 
transfer  as  well  as  those  indicated  above.  A  closed  form 
solution  for  e  was  sought  for  both  1-3  exchangers  with  one 
having  two  out  of  three  tube  passes  In  parallel  flow  and  the 
other  having  two  out  of  three  tube  passes  in  counterfiow. 

The  analytical  approach  taken,  and  demonstrated  in  this 
section,  is  for  two  out  of  three  tube  passes  in  counterflow. 
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B.  ANALYTICAL  DEVELOPMENT 

The  derivation  for  the  effectiveness,  e,  of  the  1-3 
exchanger  as  a  function  of  the  capacity  rate  ratio,  R,  and 
number  of  transfer  units,  Ntu»  depends  on  several 
assumptions. 

(1)  The  overall  coefficient  of  heat  transfer,  U,  does 
not  vary  within  the  exchanger. 

(2)  The  specific  heat  of  both  hot  side  and  cold  side 
fluids  does  not  vary. 

(3)  Each  fluid  is  thoroughly  mixed,  that  is,  the 
temperature  of  both  hot  and  cold  side  fluids  is 
uniform  over  any  cross  section. 

(4)  Steady  flow  conditions  are  maintained. 

(5)  Heat  losses  to  or  from  the  environment  are 
negligible. 

(6)  No  change  of  phase  takes  place;  all  heat 
transferred  is  sensible  heat. 

(7)  There  is  equal  heat  transfer  surface  in  each  pass. 

The  configuration  is  shown  in  Figure  3.1  where  the  three 

tube  passes  are  designated  with  subscripts  a,  b  and  c.  The 
temperature  of  the  hot  (shell  side)  fluid  is  indicated  by 
upper  case  letters.  For  the  cold  (tube  side)  fluid,  lower 
case  letters  are  used.  The  subscript  1  always  refers  to  the 
fluid  inlet  and  the  subscript  2  always  refers  to  the  fluid 
outlet . 

With  Wh  and  Cph  designating  mass  flow  (kg/sec)  and 
specific  heat  ( Joules/kg"*K)  of  hot  fluid  entering  at  T^  and 
leaving  at  T2  we  define  a  capacity  rate  for  the  hot  side 


^h  *  ^hCph 


Development  of  Effectiveness  Relationship 


.  *.  *.  v  v  r. 


In  similar  fashion  for  the  cold  side  (with  Wc  and  Cpc) 
entering  at  ti  and  leaving  at  t2>  we  have 

Cc  “  WcCpc 

We  then  obtain  an  energy  balance  for  the  entire 
exchanger 

Ch<Tl  -  T2)  -  Cc(t2  -  ti)  (1 

Over  the  right  hand  side  of  the  exchanger  (Figure  3.1) 
ch(Ti  -  T)  -  Cc(t2  -  tc  ♦  tb  -  ta)  (2 

and  a  differentiation  gives 

ChdT  »  Cc(dtc  -  dtb  +  dta)  (3] 

Across  dx,  with  a  (m2/m),  the  surface  per  running  meter 
of  length  of  pass  so  that  A  -  3aL  is  the  total  surface  in 
the  exchanger,  we  may  write  the  heat  transferred  to  the 
element  dx  in  each  cold  pass. 


Ccdta  * 

Uadx(T 

"  ta) 

(4a 

Cc^tb  * 

-Uadx(T 

-  tb) 

(4b 

Cc^c  m 

Uadx(T 

—  tc ) 

(4c 

Here  it  should  be  observed  that  due  cognizance  has  been 
taken  of  the  direction  of  the  flow  in  each  cold  fluid  pass 


with  respect  to  the  positive  sense  of  the  length  coordinate, 
a,  and  D  is  the  overall  heat  transfer  coefficient  (W/m^-'C). 
With  eqs.  (4)  in  eq.  (3) 

CfcdT  *  Ua(3T  -  ta  -  tb  -  tc)dx 


or 


dT 

Ha 


nh(3T  "  *a 


-  t. 


-  V 


(5) 


where 


is  a  sort  of  Ntu  per  unit  length  for  the  hot  side. 
Now  differentiate  eq.  (5) 


and  with  eqs.  (4)  substituted 

2 

T?  “  3nh  Hi  "  ncnh(T  "  ta  +  tb  “  tc)  (6) 


where 


n 


c 


Ua 

XT 

c 


where  again  the  resemblence  of  nc  to  N^u  can  be  noted. 


From  eq.  (2)  we  obtain 


C” 

c 


(T  -  T)  -  t„  -  t.  -  t  -  t 


and  with  eq.  (7)  put  into  eq.  (6)  we  obtain 


~n  dT 
3nh  3? 


“  n^n.  [T  + 


TT  <T1  -  T> 


where 


r?  -  3nh  §  -  -no°h  rr  f<Rc  -  1>T  +  Ti  -  Vs'  (l,> 

ax  c 


Rc  “  ^-c/^h 


is  the  capacity  rate  ratio  for  the  cold  side. 
Notice  that 


Da  %  Ua  #  Ch  .Ua  2 


ncnh  *  TT~  'XT'  *  TT~ 


(6=)  =  m 
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R  -  1  »  °h 

Rh  "  TT  “  tr 

c  c 


a  capacity  rate  ratio  for  the  hot  side.  Then,  algebraic 
adjustment  provides 


which  is  a  linear,  non-homogeneous ,  second  order 
differential  equation  with  constant  coefficients  having  a 
complementary  function 


T 

c 


a1  x  cu>: 

C,e  +  C„e 


( 


where  Ci  and  C2  are  arbitrary 


constants  and  where 


a 
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a2 


1/2 


But 
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XT 
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so  that 


where 


A  -  [9  -  41^(1  -  V}1/2  (12) 

Designate  the  particular  integral  as  Tp*  and  by  the 
method  of  undetermined  coefficients  let  Tpi  *  P  so  that  in 
eq.  (9) 


1  - 


m(- 


*h 


)  P  -  m(^  -  Tl) 


This  makes 


pi 


p  -  '  ^  -  T!l  'r^> 


so  that 


(13) 


The  general  solution  to  eq.  (9)  is  the  sum  of  eqs.  (10)  and 
(13) 


T(x) 


OjX  x  tp  —  j 

Cle  +  C2e  +  "f  -  TT 


(14) 


where  the  arbitary  constants,  C?i  and  C2  are  evaluated  from 
conditions  at  x  *  0  and  x  ■  L.  At  x  *  0,  T(x  *0)  *  T2  and 
at  x  *  L,  T(x  ■  L)  »  Ti .  When  these  are  inserted,  in  turn, 
into  eq.  (14),  one  obtains  a  pair  of  linear  algebraic 
equations  in  the  unknowns  Ci  and  C2 


where  Tpi  is  given  by  eq.  (13). 

It  is  only  a  matter  of  algebra  to  show  that 


It  is  easy  to  see  from  eq.  (1)  that 


so  that 

t2  “  ti  +  Rh(Tx  -  T2> 

Use  of  this  in  eq.  (13)  shows  that 

t,  +  Ru(T.  -  T„ )  -  R.T. 


indicating  two  alternative  forms  for  Tpi  given  by  eqs.  (13) 
and  (16). 


and  if  we  subtract  eq.  (4a)  from  eq.  (4c)  we  obtain 


dt  -  dt 
a  c 


*a  -  *c 


rdxs  _ncdx 


which  can  be  integrated  using  C3  as  the  constant  of 
integration . 


ta  “  tc  "  C3e 


"ncx 


and  at  x  *  0  where  ta  *  ti  and  tc  =  tbc 


tl  '  tbc  *  C3 


or 


tbc  “  ti  -  C3 


In  addition  at  x  =  L,  ta  =  tab  and  tc  =  t3  so  that 


tab  _  t2  C3e 


ncL 


or 


tab  -  t2  n  L 

C3  -  --nhE  -  <*ah  ' 

tl 


This  Rives  a  relationship  between  ta^  and  t>>c 


where  Nc  *  ncL  can  be  considered  as  the  total  number  of 
transfer  units  for  the  cold  side. 

Return  now  to  eq.  (18)  and  look  at  the  conditions  at 
x  -  L  where  ta  -  tfe  *  tab,  tc  «  tz  and  T  ■  Tj_ .  These 
conditions  in  eq.  (18)  give 

nh<3Tl  "  2tab  "  t2)  =  “lCle  +  a2C2e 

where  again  we  remember  that  Cj  and  C2  are  known  constants. 
Solving  for  tab 

2tab  "  "  HT  (aicieai  +  <*£2^  >  +  3T!  '  *2 
h 

and  with  this  in  eq.  (20) 

N  1  ouL  o^L 

2tb(,  *  e  °r—  (o^e  +  OgC^e  )  +  3t2  -  3T1]  +  2tj_ 
h  *  '* 

Then  with  eq.  (21)  in  eq.  (19) 

Nc 

o1C1  +  a2C2  *  3nh  f  (T2  -  tj_)  +  e  (T 1  ~  tg)l 

N  qu.L  apL 

-e  C(aLC1e  +  a^e  ^  )  (22 

Equation  (22)  confirms  Fischer's  result  [Ref.  12: 
pp.  377-383]  and  at  this  point  in  his  development  he 
branches  off  to  seek  an  expression  for  the  Logarithmic  Mean 


Temperature  Difference  Correction  Factor,  F.  The  attention 
here  is  focused  on  e  *  f(R,  N^u)  and  the  balance  of  this 
section  continues  in  this  vein. 

Look  at  the  01C1  +  0^2  term  in  eq.  (22).  Use  of  eq. 
(11)  allows  the  representation  of  aiCi  +  a2C2 

^ —  (3  +  X)C^  +  "2 —  (3  —  X)Cg  (23a 

or 
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Then  from  eqs.  (17) 
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Moreover 
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Now  let 


so  that 
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Moreover,  it  is  easy  to  see  that 
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If  eqs.  (24)  through  (28)  are  collected  and  put  into  the 
expression  of  eq.  (23b),  the  result  is 
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(231 


^2  ~  ^1 

-  rT^V'  +  z 


2(T1  “  t2)  “  (T2  "  tl)2e 


oil 


csch  zL 
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(1  -  R^e^  sinh  zL 


-)e  csch  zL  -  (--7  _  w  ■)  coth  zL]  (29a) 


and  this  could  also  be  written  as 


^2  ~  ^1 
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The  next  step  is  to  reduce  the  right  hand  side  of  eq. 
(22).  Ose  of  eq.  (11)  permits  the  representation 


ot.  L  o to L  3nh  t  3nh  An 

al°le  +  a2C2e  =  ('~T~  +  ~2“)Cie^  +  "  ~2 


h  ^h  “l1 
-  )C2e 


or  with  eqs.  (26)  inserted 


a.L  a«L  a.L  a«L  cuL  a«L 

al^le  +  a2^2e  *  ui(C^e  +  Cge  )  +  z(C^e  -  C,,e  ) 
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But  by  eqs.  (11)  and  (17) 
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or 
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Cie  -  C2e  may  also  be  simplified.  Again  using 


eqs.  (11),  (17a)  and  (17b) 


otj^L  o^L 


a,L  a0L  (T1-t0)(e  +  e  )“(T2”t1)fe  +  e 


i  o 

Cle  “  C2e 


- snc - 

(1  -  Rh)(e  -  e 


(  (  “i4 

out 
^  ) 


The  exponential  term  at  the  far  right  in  the  numerator 
is  really  quite  simple.  From  eq.  (11) 


3n,. 


and  by  eq.  (26a)  aj  +  a2  *  3nh  ■  2w.  Thus  with  the 
combination  of  exponentials  in  the  numerator  and  the 
denominator  given  by  eqs.  (27)  and  (28)  we  find  that 

oxL  a2L  2(TX  -  tgje*4*  cosh  zL  -  (T2  -  t1)2e2aL 

C.  e  -  C0e  ™  - j - 

2(1  -  \)e  sinh  zL 


or 
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Now  with  eqs.  (31)  and  (32)  put  into  eq.  (30)  we 


obtain 
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With  eqs.  (29b)  and  (33)  inserted  into  eq.  (22)  we 
obtain 
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We  wish  to  develop  an  expression  for  the  exchanger 
effectiveness  so  we  designate  the  hot  side  effectiveness  as 


(35) 


and  we  begin  by  simplifying  eq.  (34)  by  dividing  throughout 
by  (T2  -  ti)/(l  -  Rjj )  to  obtain 
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where 


W  ■  b)L 
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and 


R  -  3nh(l  -  Rh) 
We  can  then  let 


and  get  eq.  (36)  to  look  like 
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It  is  now  a  matter  of  algebra  to  solve  for 
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The  next  step  is  to  represent  $  as  a  function  of 
This  is  done  with  some  algebraic  gymnastics  as  follows: 
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But  Rh  “  (*2  “  ti)/(Ti  -  T2)  so  that 


or 


where  $  is  given  by  eq.  (38) 

The  neatness  of  the  form  of  eq.  (39)  is  deceptive 
because,  unfortunately,  it  cannot  be  used  to  determine  a 
unique  value  for  e^.  The  reason  for  this  can  be  found  in  an 
Inspection  of  eq.  (38)  which  provides  the  value  of  4>  which 
is  used  in  eq.  (39). 

Notice  in  eq.  (38)  that  Z,  W  and  Nc  are  all  functions  of 
the  product  a L.  On  the  other  hand,  w,  z  and  R  are  functions 
of  a  only.  Thus,  it  is  impossible  to  vary  a  and  L  indepen¬ 
dently  and  still  achieve  a  unique  solution. 

For  example,  suppose  a  *  50  m^/m  and  L  *  5  m  so  that 
aL  *  250.  A  value  of  $  may  be  obtained  from  eq .  (38)  using 
these  values.  However,  if  a  *  100  and  L  =  2.5  so  that  aL  is 
still  equal  to  250,  an  entirely  different  value  of  <|>  is 
obtained  because  a  ■  100  rather  than  50. 
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One  should  resist  the  temptation  to  multiply  numerator 
and  denominator  of  eq.  (38)  by  L  thereby  creating  a 
situation  where  only  Z,  W  and  Nc  appear  along  with  a  new 
R'  ■  RL.  Such  a  procedure  is  doomed  to  failure  because  in 
dealing  with  an  equation  derived  from  n  equations  in  n+1 
unknowns,  one  cannot  create  the  n+lth  equation  by 
multiplying  one  of  the  n  equations  by  a  constant.  This 
makes  the  n+lth  equation  so  obtained  linearly  dependent  on 
one  of  the  original  n  equations  and  the  entire  set  becomes 
linearly  dependent. 

This  section  represents  an  attempt  to  obtain  e  =  f(R|Ntu) 
and  the  attempt  has  not  been  successful.  It  is  now  time  to 
turn  to  the  computer  and  this  will  be  done  in  Section  IV. 


NUMERICAL  AND  COMPUTER  ANALYSIS 


With  a  closed  form  solution  for  e  as  a  function  of  R 
and  Ntu  not  attainable  as  indicated  in  Section  III,  it 
becomes  apparent  that  an  alternative  method  is  needed  to 
determine  the  effectiveness  for  the  1-3 :2C  and  1-3 :2P  heat 
exchangers.  Kern  and  Kraus  [Ref.  12:  pp.  306-360]  describe 
a  computer  code  for  a  thermal  analyzer.  This  code  (program) 
makes  use  of  node  equations  generated  by  finite  differences 
and  it  employs  a  Cholesky  LU  decomposition  scheme. 

The  Cholesky  decomposition,  as  explained  by  Stewart 
[Ref.  13:  pp.  134-144],  is  best  used  when  decomposition  in 
the  presence  of  positive  definite  matrices  is  requested. 

This  is  the  case  at  hand  where  one  tries  to  solve 
numerically  for  the  temperatures  that  lead  to  the 
effectiveness  of  the  1-3 :2C  and  1-3 :2P  heat  exchangers. 

A.  THERMAL  ANALYZER  TVSSI 

The  computer  program  employed  is  Program  TVSSI  (Appendix 
A)  which  is  an  adaptation  of  the  thermal  analyzer  program 
called  TVSS2  and  listed  by  Kern  and  Kraus  [Ref.  12:  pp. 
306-360].  The  adaptation  consisted  of  changing  the  program 
so  that  it  could  perform  the  computations  in  the  SI  system 
of  units  and  be  receptive  to  the  use  of  a  specially  created 
input  file.  Also  it  should  be  noted  that  TVSS2  was  written 


to  be  used  in  conjunction  with  the  Honeywell  H-1800  computer 
system.  Therefore,  it  had  to  be  modified  to  run  on  the 
Naval  Postgraduate  School's  IBM  3033  AP  system. 

The  program  itself  is  a  non-linear  equation  solver  that 
determines  the  temperatures  at  a  prescribed  number  of  node¬ 
points  or  nodes  from  a  set  of  node  equations  in  almost  any 
framework  (i.e.,  network  analysis,  field  plotting,  or  fluid 
flow  distribution).  It  has  certain  features  that  make  it 
primarily  an  equation  solver  for  thermal  analysis.  These 
features  include: 

1.  an  ability  to  linearize  radiation  terms. 

2.  an  ability  to  allow  any  of  the  coefficients  in  the 
node  equations  to  vary  with  temperature. 

3.  an  ability  to  provide  constant  heat  input  and  heat 
input  as  a  function  of  temperature  at  any  node. 

4.  an  ability  to  consider  other  modes  of  heat  transfer 
that  are  non-linear  such  as  boiling  and  natural 
convection . 

As  stated  earlier,  the  program  utilizes  the  Cholesky 
decomposition  scheme  and ,  because  of  the  linearization  of 
the  radiation  terms  (a  feature  of  the  program  that  is  used 
even  though  radiation  does  not  appear  in  this  e  -  Ntu 
study),  the  program  is  iterative. 

Cholesky' s  decomposition  consists  of  finding  a  lower 
triangular  matrix  [L]  which  is  capable  of  reducing  the 
original  system  of  equations. 


or 


mm  -  fQ]  -  [o]  (4D 

to  the  unit  triangular  form 

UHT]  -  fB]  *  JO]  (42) 

so  that  the  sought  after  elements  of  the  column  vector  fT] 
can  be  obtained  by  backward  substitution. 

Suppose,  for  example,  that  [K]  is  3x3  and  assume  that 
the  system  [K][T]  *  [Q]  has  been  reduced  to  the  form 
fA][T]  -  [B]  ■  [0].  In  this  event  a  premultiplication  by 
[L]  will  return  the  system  to  its  original  form,  that  is 

[L]([A][T]  -  [B] )  *  [ K ] t T ]  -  [Q]  =  [0] 

This  implies  that 

[L][A]  *  [K]  (43) 

and 

f L] fB]  »  rO]  (44) 

These  equations  allow  the  determination  of  fL],  f A ] ,  and 
fB]  in  a  very  simple  manner  and  the  matrices  are  uniquely 
determined  because  fK]  and  fQ]  are  known,  or,  at  least  are 
known  after  each  iteration  because  the  elements  of  ]K]  are 
linearized.  For  a  3x3  system 
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one  may  obtain  the  following 
and  [ B ] . 
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the  elements  of  [L],  [A], 


kil  *  (l)lil  +  (0)li2  +  (0)li3  ■  lil  (45) 

which  shows  that  the  first  column  of  f L ]  is  identical  to  the 
first  column  of  [K]. 


klj  =  ( 1 1 1 ) al j  +  (°)a2j  +  (°)a3j  “  1Haij  =  kllai.j  (46> 
which  shows  that  the  first  row  of  fA]  is  equal  to  the  first 


row  of  f K ]  divided  by  fku]  and 
k22  =  121a12  +  122 ( 1 )  * 
k23  *  !21a13  +  !22a23  * 
k32  “  !31a12  +  132<1>  * 
Q2  =  l21bl  +  ^22b2 
In  the  foregoing  manner  the 


then 

122  *  k22  -  l21a12 
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elements  of  FL],  [A]  and  fRl 


are  obtained  successively  in  terms  of  previously  determined 
elements  in  a  progression  that  goes  horizontally  from  I22 
on.  Thus  the  general  relationship  are  seen  to  be 
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Xij  “  kij  I  xirkrJ 
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(47) 


with 


lil  *  kil 


(48) 


and 


(49) 


with 


Moreover,  it  is  observed  that  if  TK1  is  symmetrical  which  it 
must  be  in  our  coupled  set  of  equations  (kj.j  “  kji)  then 

lj.  j  *  a-jilii  ( i  ,  j  *  1,2,3,...  n— 1 ;  i*j)  (51) 

The  modification  for  computation  in  the  SI  system  was 
quick  and  simple.  It  involved  changing  a  numeral  in  two 
places  (460  to  273)  and  some  format  statements  (*F  to  °C  and 
Btu/hr  to  Watts). 

The  conversion  of  TVSS2  to  TVSSI  for  running  on  the 
IBM  3033  AP  system  (which  uses  the  FORTVS  compiler)  required 
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a  modification  to  the  Fortran  program  language  used  in  TVSS2 
in  order  to  compile  the  FORTVS  (basically  international 
Fortran  77)  system  used  on  the  IBM  3033. 

B.  INITIAL  MODELING 

Initially,  the  model  was  designed  to  find  the  tempera¬ 
tures  that  would  allow  computation  of  a  single  effectiveness 
value  after  a  detailed  set  of  capacity  rate,  coefficients 
and  surface  data  were  entered  into  the  program.  From  the 
scope  of  the  problem  it  was  realized  that  multiple  runs  of 
the  thermal  analyzer  program  (TVSSI)  would  be  needed.  It 
therefore  became  necessary  to  develop  a  program  that  given 
Ch»  Cc,  U,  A,  Ti ,  and  ti ,  an  input  file  would  be  created  for 
use  by  the  modified  version  of  the  thermal  analyzer 
(TVSSI). 

The  first  step  taken  was  to  develop  a  program  to  create 
an  input  file  for  TVSSI  that  would  yield  the  effectiveness 
for  a  1-4  heat  exchanger  which  could  be  compared  to  the 
existing  analytical  solution  for  the  effectiveness  of  a  1-4 
exchanger.  With  this  accomplished  and  confidence 
established,  similar  programs  for  the  1-3:2C  and  1-3:2P 
exchangers  could  be  developed.  This  program  was  called 
NTU14  (See  Appendix  B)  and  the  following  parameters  were 
used  for  all  runs. 

1.  250  nodes  were  used. 

2.  The  initial  temperature  for  the  computer  to  begin  the 
iterative  process  was  set  at  200°C. 
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3.  An  eventual  accuracy  of  .05  between  the  final  and  next 
to  last  iterations  was  used. 

4.  A  radiation  coefficient  convergence  factor  of  0.66667 
between  iterations  was  used. 

5.  The  maximum  number  of  iterations  that  the  computer 
was  allowed  to  perform  was  set  at  12. 

6.  A  damping  factor  of  .8  was  set  as  an  initial  damping 
based  on  the  number  of  non-linear  terms  in  all  of  the 
node  equations. 

When  the  values  of  Ch,  Cc,  U,  A,  and  ti  are  set,  N^u 
and  R  are  then  compiled  and  an  input  file  for  TVSSI  was 
generated.  In  this  file  all  node  equations  and  internode 
conductance  values  were  determined.  This  program  determined 
and  specified  the  nodes  that  interact  with  each  other  and 
the  methods  by  which  the  interaction  takes  place  such  as 
conduction,  forced  convection,  and  fluid  flow. 

The  program,  NTU14 ,  makes  use  of  the  fact  that  each  term 
in  a  node  equation  shows  three  things.  The  first  is  the 
node  that  is  coupled  for  heat  flow  with  the  node  in 
question.  The  second  is  the  method  of  heat  flow  between  the 
nodes.  In  this  case  forced  convection  and  fluid  flow  are 
used.  Finally,  the  node  equation  shows  the  magnitude  of  the 
internode  heat  flow.  Here  all  the  pieces  of  information  are 
collected  and  presented  for  use  by  TVSSI  as  an  input  file 
with  all  items  in  the  proper  format. 

A  comparison  of  the  effectiveness  for  the  1-4  exchanger 
developed  by  the  computer  to  that  of  using  the  closed  form 
analytical  solution  for  effectiveness  developed  by  Kraus 


and  Kern  [Ref.  8]  as  shown  by  equation  (52) 


1  +  R  +  i[l  +  4R2]1/2coth(  foL1  *  ^-1)  +  4  tanh 


was  then  undertaken.  The  results  of  this  comparison  showed 
that  over  the  entire  range  of  R  from  .01  to  1.0  for  varying 
values  of  N-tu  from  0  to  3.25  less  than  a  0.5%  difference  was 
ever  realized.  A  small  sample  of  these  results  are  provided 
in  Table  1.  The  conclusion  to  be  drawn  here,  is  that  the 
methodology  used  to  develop  the  computer  program  NTU14  for 
input  to  TVSSI  for  finding  effectiveness  was  sound  and  could 
then  be  used  in  the  development  of  the  1-3 :2C  and  1-3 :2P 
exchanger  methodology. 


C.  DEVELOPED  MODELS  FOR  1-3 :2C  AND  1-3 :2P  HEAT  EXCHANGERS 
The  same  technique  used  in  developing  the  program  NTU14 
was  used  to  generate  computer  programs  NTU32C  and  NTU32P . 
These  are  listed  in  Appendices  C  and  D.  The  departure  for 
each  of  these  programs  from  the  NTU14  program  is  in  the 
number  of  nodes;  they  are  based  on  200  node  models  as  shown 
in  Figures  4.1  and  4.2.  An  example  of  the  output  file 
generated  from  one  of  these  programs  is  found  at  Appendix  E. 
It  is  these  values  shown  in  Appendix  E  that  are  used  by  the 
thermal  analyzer  to  determine  the  temperatures  T2  and  t2 
for  the  specific  set  of  given  initial  parameters  ,  Cc ,  R, 
A,  Ti  ,  and  ti. 


TABLE  1 

COMPUTER  TO  ANALYTICAL  COMPARISON 
FOR  R  =  .5 


D.  SCOPE  OF  COMPUTER  ANALYSIS 

At  this  point  it  is  possible  to  let  the  computer  solve 
for  temperature  values  that  yield  a  value  for  effectiveness 
based  on  a  specific  set  of  initial  parameters.  However, 
it  must  be  realized  that  many  computer  runs  are  required  to 
generate  enough  data  to  ensure  confidence  in  the  results 
which  cover  a  wide  range  of  capacity  rate  ratios  and  Ntu 
values. 

To  efficiently  expedite  the  computer  task,  the  Multiple 
Virtual  System  (MVS)  with  Job  Entry  Subsystem  and  Networking 
(JES3)  was  utilized.  The  MVS  coupled  with  JES3  is  more 
commonly  referred  to  as  batch  processing.  Based  on  trial 
and  error,  it  was  determined  that  in  order  to  build  a  solid 
data  base,  eleven  different  values  for  effectiveness  were 
needed  to  best  represent  a  particular  value  of  R.  This 
was  required  over  a  range  of  R  from  R  *  0.1  to  1.0  in 
increments  of  .01.  In  all,  200  curves  for  both  the  1-3 :2C 
and  1-3 :2P  exchangers  were  needed.  This  means  that  2,200 
unique  values  of  effectiveness  needed  to  be  found,  plus 
100  values  of  effectiveness  to  be  used  for  comparison  with 
the  1-4  exchangers. 

To  complete  this  task,  TVSSI  was  slightly  modified  in 
accordance  with  the  appropriate  guidelines  of  the  ,1ob 
control  language  (JCL)  needed  to  run  on  the  batch  processing 
system.  These  modifications  are  few  and  were  needed  only 
at  the  beginning  and  end  of  TVSSI.  The  modified  version 


of  TVSSI  has  been  called  TVCOUNT  with  changes  shown  in 
Appendix  F.  It  was  TVCOUNT  that  was  then  used  to  activate 
TVSSI . 

It  also  became  necessary  to  modify  the  three  input  file 
programs  NTU14,  NTU32P  and  NTU32C  such  that  they  needed  to 
be  compiled  only  once.  They  were  then  loaded  in  a  library 
file  to  be  used  when  called  by  another  program.  New 
programs  utilizing  the  batch  system  were  written  that  could 
easily  be  loaded  with  the  appropriate  input  data  for  a 
specific  R  value.  These,  which  are  referred  to  as  "sister 
programs,"  were  used  to  go  from  the  library  file  to  TVSSI 
and  cause  TVSSI  to  be  executed  eleven  times  under  TVCOUNT 
covering  the  desired  range  of  N-tu  for  a  specific  R  value. 

The  revised  input  files  called  NTU14C,  NTU32CC ,  NTU32PC  and 
their  associated  "sister  execution  programs,"  NTU14L , 

NTU32CL ,  NTU32PL ,  are  found  in  Appendices  G  through  L. 

The  overall  system  flow  chart  of  how  all  of  the  fore¬ 
going  is  accomplished  is  found  in  Figures  4.3,  4.4  and  4.5. 
It  is  noted  from  these  figures  that  TVSSI  is  referred  to  as 
TVSSIA  through  TVSSIV.  These  are  the  same  programs  as  TVSSI 
but  for  bookkeeping  purposes  by  the  computer  they  are 
labeled  A  through  V. 

E.  COMPUTER  RESULTS 

Upon  completion  of  all  data  collection  from  the  computer 
output,  plots  of  effectiveness  vs.  N-tu  for  the  whole  range 
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of  R  were  plotted  and  are  shown  In  Appendices  II  and  N. 
Thirty-three  different  plotting  programs  utilizing  the 
"Display  Integrated  Software  System  and  Plotting  Language 
(DISSPLA),"  were  written  to  graph  the  data  obtained. 

An  example  of  one  of  these  programs  is  provided  in 
Appendix  0. 

Examination  of  the  graphs  in  Figures  4.6,  4.7,  and  4.8 
shows  that  the  1-3 :2C  exchangers  outperform  the  1-3 :2P,  1-2 
and  1-4  exchangers.  This  is  true  in  all  cases,  and,  because 
of  this,  only  a  sample  of  the  data  was  chosen  to  be  shown  in 
these  figures.  Furthermore,  it  is  noted  that  at  higher  N^u 
values,  the  effectiveness  of  the  1-3 :2C  exchanger  is  better 
than  all  of  the  others  considered,  while  at  higher  capacity 
rate  ratios,  the  effectiveness  of  1-3 :2C  exchanger  even 
begins  to  outperform  the  others  at  lower  Ntu  values.  This 
increase  in  performance  is  easily  understood  because  it  has 
been  proven  by  Kern  [Ref.  1:  pp.  139-137]  and  others  that 
greater  temperature  differences  result  when  process  streams 
are  in  counterflow  than  parallel  flow.  Thus,  when  there  is 
a  combination  of  the  two  phenomena  (counterflow  and  parallel 
flow)  occurring,  then  it  becomes  apparent  that  the  extra 
counterflow  pass  must  increase  the  exchanger's  overall 
performance . 

As  shown  in  the  graphs  in  Appendices  M  and  N,  effective¬ 
ness  increases  as  both  R  and  N^u  increase.  These  curves  can 
be  used  to  give  a  graphical  approximation  of  effectiveness 
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if  that  is  all  that  is  required.  Prom  the  empirical  data 
used  to  develop  these  curves,  further  investigation  into  the 
development  of  equations  for  these  curves  which  may  be  used 
to  determine  an  exact  value  of  effectiveness  when  R  and  N^u 
are  known,  is  undertaken  in  Section  V. 
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11  sets  of  initial  parameters 
given  (Ntu,  Cj,.  Cc,  U,  A,  Tx  and  tj) 
make  up  the  input  data  to  be 
solved  for 


Figure  4.5 


Computer  Systems  Flow  Chart  for  1-3: 2P  Exchanger 
Model 
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EFFECTIVENESS  VS.  NTU  FOR  R=1.0 


Figure  4.8  Comparison  of  Analytical  (1-2  and  1-4).  to  Computer 
Results  (1-3:2C  and  1-3:2P)  at  R  =  1.0 
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POLYNOMIAL  REGRESSION 


A.  DEVELOPMENT  OF  POLYNOMIAL  EQUATIONS 

The  empirical  data  obtained  for  the  two  1-3  heat 
exchangers  was  designed  to  cover  an  extensive  range  of  R 
values  varying  from  0.01  to  1.0  in  increments  of  0.01.  As 
discussed  earlier  in  Section  IV,  the  data  obtained  at  a 
specific  value  of  R  is  the  computer  evaluated  result  of  the 
effectiveness,  for  an  associated  N-tu  value.  With  this 
accomplished,  it  then  becomes  possible  to  graph  separate 
curves  for  each  of  the  different  R  values  as  shown  in 
Appendices  M  and  N.  Through  a  polynomial  regression 
technique,  as  discussed  in  this  section,  it  is  also  possible 
to  develop  implicit  equations  for  the  curves  with  e  =  f(Ntu, 
R).  It  is  also  apparent  from  an  inspection  of  the  graphical 
representation  of  the  empirical  data  in  Appendices  M  and  N, 
that  the  curves  conform  to  a  high  degree  polynomial. 

However,  further  analytical  investigation  is  needed  to 
ascertain  the  exact  order  of  the  polynomial  terms.  This 
investigation  will  not  only  lead  to  the  order  of  the 
polynomial,  but  to  the  specific  equation  for  each  curve. 

By  use  of  polynomial  regression,  the  least-squares 
method  can  be  readily  extended  to  best  fit  the  data  to  the 
mth-degree  for  the  polynomial 

y  *  Aq  +  Aix  +  A2X2  +  ...  Amxm  (53) 


with  the  error  defined  by 


ei  ■  Di  -  yi  »  Di  -  Ao  -  Aix  -  A2*i2  -  ...  -  AmXim 

where  represents  the  empirical  data  value  corresponding 
to  xi,  xi  being  free  of  error. 

The  objective  is  to  minimize  the  sum  of  the  squares  of 
the  residuals,  Sr, 


m  n  m 

z  e,  -  z  (D, 
i=l  1  1 


Aq  -  Aixi  "  A2x2 


Amxm)  (54) 


Because  at  a  minimum,  the  partial  derivatives  3Sr/3Ao, 
3Sr/3Ai  ...  3Sr/3Am  vanish,  after  taking  the  derivative  of 
Sr  with  respect  to  each  of  the  coefficients  of  the  poly¬ 
nomial  ,  it  can  be  seen  that 


3S 
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A  mN 
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Then  by  dividing  by  -2  and  rearranging  we  obtain 


AoM  +  Ai  I  Xi  +  Ai*i2  +  ...  +  Am  E  *  E  Di 

Aq  Z  xi  +  Ai  E  xi2  +  A2  E  Xi3  +  ...  +  Am  E  Xira+1  *  E  x^ 

A0  E  Xi2  +  Aj.  E  Xi3  +  A2  E  Xi4  +  ...  +  Am  E  xi®+2  =«  E  xi^ 

•  •  •  • 

•  •  * 

•  •  •  *  * 

A0  E  Xim  +  Ax  E  Xi^+1  +  A2  E  Xi^+2  +  ...  +  Am  E  Xi2m  -  E  Xi^Dj 

where  all  suramations  are  from  i*l  through  n.  All  of  the 
foregoing  m+1  equations  are  linear  and  have  m+1  unknowns: 

Ao,  Ai,  A2,  ...  Am.  The  coefficients  of  the  unknowns  can  be 
calculated  directly  from  the  observed  data.  Thus,  the  pro¬ 
blem  of  determining  a  least-squares  polynomial  of  degree  m 
is  equivalent  to  solving  a  system  of  m+1  simultaneous  linear 
equations.  Putting  the  equations  in  matrix  form  yields 


N  E  Xi  E  xi2  E  Xi3  .  .  .  E  Xim 

E  Di 

E  Xi  E  Xi2  E  Xi3  E  Xi4  •  •  ♦  E  Xin+1 

E  XiDi 

E  Xi2  E  Xi3  E  Xi4  E  Xi5  .  .  .  E  Xim+2 

#  •  •  •  • 

[A]  = 

E  Xi2Di 

• 

•  •  •  •  • 

•  t  •  •  • 

E  Xim  E  Xim+1  E  Xim+2  E  Xim+3  .  .  .  E  xt2m 

' 

• 

E  XimDi 

_  - 

[Ref.  15:  pp.  302-309  and  Ref.  16:  468-4741. 

From  this  point  on,  one  finds  that  it  is  best  to  use  a 
computer  to  assist  in  solving  the  simultaneous  equations  and 
this  will  also  help  alleviate  any  ill-conditioning  that  may 
otherwise  occur.  An  existing  "curvefit"  program  available 
through  NON-IMSL  (Ref.  16]  and  found  in  Appendix  P  was 
used  although  some  modifications  were  made  to  the  original 
program  to  best  accommodate  the  goals  of  this  work. 

To  determine  the  order  of  polynomial  that  should 
eventually  be  used,  one  increases  the  degree  of  the 
approximating  polynomial  as  long  as  there  is  a  statistically 
significant  decrease  in  the  variance  o^,  which  is  computed 
by 
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In  otherwords,  the  selection  of  the  optimum  degree  poly¬ 
nomial  is  contingent  upon  a  decreasing  variance  and  once  the 
variance  begins  to  increase,  the  degree  of  the  polynomial 
becomes  too  high.  For  all  cases,  it  was  found  that  the 
e  -  N^u  developed  curves  are  of  the  5th  order. 

As  shown  in  Figures  5.1  and  5.2  the  computed  values  of 
effectiveness  vs.  N^u  for  R  *  0.1,  0.5  and  1.0  for  both 
flow  arrangements,  (1-3:2P)  and  (1-3:2C),  have  been  graphed 
and  fitted  by  a  5th  degree  polynomial.  Because  all  computed 
values  for  effectiveness  follow  a  predictable  trend, 
only  a  sample  of  the  data  covering  the  whole  range  of 


values  of  R  have  been  shown.  It  Is  clear  that  the  graphic 
interpretation  strongly  backs  what  Is  known  analytically 
from  the  polynomial  regression  technique.  Where  the 
relationship  for  c  »  f(Ntu.  is  found  explicitly  from 

e  *  AsNtu5  +  A4Ntu4  +  A3Ntu3  +  A2Ntu2  +  AlNtu  +  A0  (56) 

while  the  corresponding  coefficients  A5,  A4,  A3,  A2,  Ai ,  and 
A0  relating  to  a  specific  value  of  R  are  found  in  Tables  2 
and  3  for  the  (1-3 :2P)  and  (1-3 :2C)  configurations.  An 
example  of  how  to  use  this  equation  in  a  heat  exchanger 
problem  now  follows. 

B.  NUMERICAL  EXAMPLE 

Consider  a  heat  exchanger  containing  400  m2  of  (A  * 

400  m2)  of  surface  and  operating  with  an  overall  heat 
transfer  coefficient  of  80  W/m2*C.  Cold  fluid  at  a  capacity 
rate  of  10,000  W/#C  enters  the  exchanger  at  60#C.  Hot  fluid 
at  a  capacity  rate  of  20,000  W/*C  enters  the  exchanger  at 
200°C. 

1.  Find 

The  effectiveness  ( e)  and  compute  the  hot  and  cold 
fluid  outlet  temperatures  for  the  (1-3:2C)  shell  tube  pass 
arrangement . 

2 .  Assumptions 

1)  Negligible  heat  loss  to  surroundings  and  kinetic  and 
potential  energy  changes. 


2)  Constant  thermal  and  fluid  properties  for  both 
fluids. 

3.  Analysis 

Here  Cc  -  10,000  W/«C  and  Ch  -  20,000  W/*C  this 
makes  R  ■  Cmin/Cmax  “  ^c/^h  m  (T l  “  T2)/(t2  ~  ti) 

10,000/20,000  -  0.5. 

and  Ntu  -  UA/Cc 

-  80  (400)/ (10 ,000)  -  3.2 

First,  go  to  Table  2  (page  81)  for  the  (1-3:2C)  arrangement 
with  R  -  0.5  and  find  the  coefficients 

A0  -  0.1294  x  IQ-2 

Aj  *  0.98120 

A2  =  -0.66161 

A3  -  0.27938 

A4  »  -0.66456  x  10" 1 

A5  «  0.66069  x  10-2 

Then  apply  equation  (56)  for  N^u  -  3.2 

e  *  A5Ntu5  +  A4Ntu4  +  A3Ntu3  +  A2Ntu2  +  AlNtu  +  A0  (56) 


Because  Cc  < 


and  with  Ti  -  ti  -  200  -  60  *  140 *C 

t2  -  ti  -  e(Ti  -  ti) 

-  0.769  (140) 

«  107. 7*C 

Finally,  the  outlet  cold  fluid  temperature  t2  is 

t2  *  107.7  +  ti 

-  107.7  +  60 

-  167. 7*C 

and  the  fluid  temperature,  T2  is  easily  found 


T2  »  Ti  -  0.5  (t2  -  ti) 

*  200  -  0.5  (t2  _  ti) 

*  146. 1°C 

4.  Observations 

The  primary  observation  made  here  is  that  by  using 
the  5th  order  polynomial  equation  (56)  with  the  appropriate 
coefficients  found  in  Table  2  or  3,  an  accurate  value  for 


effectiveness  can  be  found  thus  allowing  one  to  solve  for 
many  more  unknown  values  or  prarameters  of  the  heat 
exchanger  (i.e.,  hot  and  cold  outlet  temperatures). 

The  other  observation  that  is  to  be  made  is  that 
when  comparing  the  value  for  effectiveness  computed  here 
against  the  value  for  a  1-2  or  1-4  exchanger  (0.745  and 
0.740  respectively)  under  the  same  conditions,  one  finds 
that  there  is  a  significant  difference  in  exchanger 
performance  as  a  function  of  odd  or  even  tube  passes  and 
that  the  1-3 :2C  arrangement  has  a  higher  effectiveness  than 
either  the  1-2  or  1-4  arrangement.  From  inspection  of  the 
curves  for  the  1-3 :2P  exchanger  at  Figure  N-l  or  N-6  with 
R  =  0.5  and  Ntu  =  3.2  an  approximate  value  of  e  »  .715  is 
obtained.  It  is  clear  that  this  is  also  less  than  that  of 
1-3 :2C  arrangement.  Therefore,  it  is  evident  that  the 
1-3 :2C  exchanger  out-performs  not  only  the  1-2  and  1-4 
arrangement  but  also  its  counterpart  the  1-3 :2P  exchanger  by 
3.1%,  3.8%  and  7.0%  respectfully. 
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TABLE  3  (cont’d) 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Effectiveness  values,  though  not  analytically  derived, 
can  be  determined  for  the  1-3 :2C  and  1-3 :2P  heat  exchangers 
by  utilizing  a  5th  order  polynomial  approximation. 

Sufficient  data  now  exists  to  cover  a  complete  range  of 
capacity  rate  ratios  for  values  of  Ntu  from  0.0  to  3.25. 

From  the  knowledge  of  a  particular  three  tube  pass  heat 
exchanger  arrangement  and  dimensions,  fluid  flow  rates  and 
temperature  extremes,  Ntu  an<*  R  values  may  be  computed. 

Then  the  effectiveness  may  be  determined  using  the 
appropriate  coefficients  from  the  tables  provided  herein. 
After  the  effectiveness  is  obtained,  there  are  two  choices. 

1)  From  a  knowledge  of  qmax  (see  Section  III),  it  is  a 
simple  matter  to  determine  the  actual  heat  transfer 
rate  from  the  postulation  that  q  =  eqmax. 

2)  Both  fluid  outlet  temperatures  may  be  computed  since 
q  =  wcCpj1(t2  ~  ti). 

It  is  also  apparent  from  the  work  done  here  that  the 
1-3 :2C  exchanger  effectiveness  outperforms  that  of  the  1-2, 
1-4  and  its  counter  part,  the  1-3 :2P  exchanger,  as  Ntu 
increases.  This  is  true  for  any  and  all  values  of  R. 
Therefore,  because  it  is  possible  to  determine  the  effec¬ 
tiveness  of  a  1-3  exchanger  which  has  a  higher  effectiveness 
(the  1-3 :2C  arrangement)  than  that  of  the  1-2,  1-4  and 


1-3 :2P  exchanger  particularly  at  high  N^u  levels,  the  1-3 
exchanger  can  now  be  given  full  consideration  in  heat 
exchanger  design.  This  would  be  especially  helpful  where 
three  tube  passes  could  alleviate  a  configuration  problem. 


B .  RECOMMENDAT I ONS 

The  following  recommendations  are  provided  for  possible 
follow-on  projects  of  a  similar  nature. 

•  Continue  study  of  the  e  -  N-tu  method  for  the  1-5  heat 
exchangers. 

•  Investigate  the  possibility  of  using  a  linear 
approximation  from  the  1-3  data  to  find  a  periodicity 
of  R  values  at  which  to  develop  a  1-5  data  base. 

•  Develop  interactive  software  to  be  used  on  a 
microcomputer  that  contains  both  the  1-3  and  1-5 
polynomial  approximation  coefficients.  The  procedure  of 
entering  the  values  for  N^u,  R,  and  the  type  of  heat 
exchanger  when  asked  to  do  so  in  a  menu-driven  fashion 
so  that  effectiveness  values  can  be  readily  obtained 

is  self-evident. 

•  After  results  of  1-5  exchanger  analyses  are  complete, 
investigate  the  effectiveness  for  both  the  1-3  and  1-5 
heat  exchangers  experimentally  to  confirm  or  refute 
the  theoretical  work  done  here. 
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APPENDIX  B 

NTU14  COMPUTER  GENERATED  INPUT  ANALYZER  PROGRAM 
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APPENDIX  D 

NTU32P  COMPUTER  GENERATED  INPUT  ANALYZER  PROGRAM 
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APPENDIX  E 

SAMPLE  OUTPUT  FROM  NTU14  COMPUTER  INPUT  ANALYZER  PROGRAM 
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APPENDIX  F 

MODIFIED  SECTIONS  OF  THERMAL  ANALYZER  TO  RUN  ON  BATCH  SYSTEM 
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APPENDIX  K 

MODIFIED  NTU32P  PROGRAM  NTU32PC  USED  TO  RUN  ON  BATCH  SYSTEM 
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UNNORMALIZE"  THE  RETURNED  COEFFICIENTS. 


REAL  FUNCTION  AMIN (X , NITEMS ) 

PURPOSE:  FIND  THE  MINIMUM  VALUE  IN  AN  ARRAY. 


INTEGER  NITEMS 


LIST  OF  REFERENCES 
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